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Abstract. We extend the formalism of the relativistic thermal and kinematical Sunyaev-Zeldovich effects to the 
observer's system (the Solar System) moving with a velocity (3s = vs/c with respect to the cosmic microwave 
background radiation. The present formulation makes full use of the Lorentz covariance properties of the problem 
and gives solutions with high precision. We confirm the results recently obtained by Chluba, Huetsi, and Sunyaev 
in the lowest order of the observer's velocity /3s ■ We give a more general analytic expression for the thermal and 
kinematical Sunyaev-Zeldovich effects corresponding to the observer's system (the Solar System) with the power 
series expansion approximation in terms of 8 e = ksTe/mc 2 , where T e and m are the electron temperature and 
the electron mass, respectively. 
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1. Introduction 

Compton scattering of the cosmic microwave background 
radiation (CMBR) by hot intracluster gas — the Sunyaev- 
Zeldovich effect (Zeldovich & Sunyaev 1969; Sunyaev 
& Zeldovich 1972, 1980a,b, 1981) — provides a useful 
method for studies of cosmology (see recent excellent 
reviews: Birkinshaw 1999; Carlstrom, Holder, & Reese 
2002). The original Sunyaev-Zeldovich formula has been 
derived from a kinetic equation for the photon distribution 
function taking the Compton scattering by electrons into 
account: the Kompaneets equation (Kompaneets 1957; 
Weymann 1965). The original Kompaneets equation has 
been derived with a nonrelativistic approximation for the 
electron; however, recent X-ray observations have revealed 
the existence of many high-temperature galaxy clusters 
(Tucker et al. 1998; Markevitch 1998; Allen et al. 2001; 
Schmidt et al. 2001; Allen et al. 2002). Some galaxy clus- 
ters have been found to possess the electron temperature 
k B T e ~20keV. Wright (1979) and Rephaeli and his collab- 
orator (Rephaeli 1995; Rephaeli & Yankovitch 1997) have 
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done pioneering work taking the relativistic corrections to 
the Sunyaev-Zeldovich effect into account for clusters of 
galaxies. 

In recent years remarkable progress has been made 
in theoretical studies of the relativistic corrections to 
the Sunyaev-Zeldovich effects for clusters of galaxies. 
Stebbins (1997) generalized the Kompaneets equation. 
Itoh, Kohyama & Nozawa (1998) have adopted a relativis- 
tically covariant formalism to describe the Compton scat- 
tering process (Berestetskii, Lifshitz, & Pitaevskii 1982; 
Buchler & Yueh 1976), thereby obtaining higher-order rel- 
ativistic corrections to the thermal Sunyaev-Zeldovich ef- 
fect in the form of the Fokker-Planck expansion. In their 
derivation, the scheme to conserve the photon number 
at every stage of the expansion proposed by Challinor & 
Lasenby (1998) played an essential role. 

Nozawa, Itoh & Kohyama (1998b) have extended their 
method to the case where the galaxy cluster is moving 
with a peculiar velocity with respect to CMBR. They have 
thereby obtained the relativistic corrections to the kine- 
matical Sunyaev-Zeldovich effect. Challinor & Lasenby 
(1999b) then confirmed the correctness of the result ob- 
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tained by Nozawa et al. (1998b). Sazonov & Sunyaev 
(1998a, b) calculated the kinematical Sunyaev-Zeldovich 
effect by a different method. 

Itoh, Nozawa & Kohyama (2000) have also applied 
their method to calculation of the relativistic corrections 
to the polarization Sunyaev-Zeldovich effect (Sunyaev & 
Zeldovich 1980b, 1981). They thereby confirmed the re- 
sult of Challinor, Ford & Lasenby (1999a) which was ob- 
tained with a completely different method. Recent work 
on the polarization Sunyaev-Zeldovich effect include Audit 
& Simons (1999), Hansen & Lilje (1999), and Sazonov & 
Sunyaev (1999). 

Sazonov & Sunyaev (1998a, b) reported the results of 
the Monte Carlo calculations on the relativistic corrections 
to the Sunyaev-Zeldovich effect. In Sazonov & Sunyaev 
(1998b), a numerical table that summarizes the results of 
the Monte Carlo calculations was presented. Nozawa et al. 

(2000) presented an accurate analyic fitting formula that 
has a high accuracy for the ranges 0.00 < 9 e < 0.05 and 
< X < 20, where 9 e = k B T e /m e c 2 , X = hu/k B T . 
Another fitting formula with a still higher precision that 
is valid for the more limited ranges 0.00 < 6 e < 0.035, 
< X < 15 was developed by Itoh et al. (2000b). 
Relativistic corrections to the double scattering effect on 
the Sunyaev-Zeldovich effect was calculated by Itoh et al. 

(2001) . Dolgov et al. (2001) carried out an independent nu- 
merical calculation, and their result shows excellent agree- 
ment with that of Itoh et al. (2001) for the case of small 
optical depth. 

Very recently Kitayama et al. (2004) measured the 
Sunyaev-Zeldovich effect in the galaxy cluster RX J1347- 
1145. They have discovered that the south-east excess 
component of this galaxy cluster has a temperature 
ksTex = 28.5 ± 7.3keV. They attribute the high temper- 
ature of this excess component to a recent major merger 
as discussed by Sarazin (2003). Therefore it is important 
to present accurate numerical data for k B T e > 25keV. 
Itoh & Nozawa (2004) have accordingly presented an ac- 
curate numerical table for the relativistic corrections to 
the Sunyaev-Zeldovich effect for clusters of galaxies in the 
ranges 0.002 < 9 e < 0.100. For analyses of the galaxy 
clusters with extremely high temperatures, the results of 
the calculation of the relativistic thermal bremsstrahlung 
Gaunt factor (Nozawa, Itoh & Kohyama 1998a) and their 
accurate analytic fitting formulae (Itoh et al. 2000b) will 
be useful. The nonrelativistic electron-electron thermal 
bremsstrahlung Gaunt factor has been also calculated by 
Itoh, Kawana & Nozawa (2002a). 

It now appears that all the necessary theoretical tools 
are in place for accurate analysis of the observational 
data of the Sunyaev-Zeldovich effect for galaxy clusters. 
Recently, however, Chluba, Huctsi, & Sunyaev (2005) 
point out the importance of the influence of the Solar 
System's motion with respect to the CMBR rest frame. 
Assuming that the CMBR dipolc is fully motion-induced, 
we deduce that the Solar System is moving with a ve- 
locity Ps = vs/c — 1.241 xl0~ 3 towards the direction 
{£,b) = (264.14° ± 0.15°, 48.26° ± 0.15°) (Smoot et al. 



1977; Fixscn et al. 1996; Fixsen & Mather 2002). Chluba 
et al. (2005) have calculated corrections to the Sunyaev- 
Zeldovich effect arising from the motion of the Solar 
System. In this paper we wish to calculate those correc- 
tions to the Sunyaev-Zeldovich effect which are due to the 
motion of the Solar System in a more general way. As a 
matter of fact, it turns out that one can make full use 
of the Lorentz covariance properties of the problem and 
obtain solutions with high precision by extending the re- 
sults obtained by Itoh, Kohyama & Nozawa (1998) and 
Nozawa, Itoh & Kohyama (1998b). 

The present paper is organized as follows. In § 2, the 
general formalism is presented in deriving the thermal and 
kinematical Sunyaev-Zeldovich effects by taking the mo- 
tion of the Solar System into account. With the power se- 
ries expansion approximation, an analytic expression in- 
cluding the relativistic effects is derived for the thermal 
and kinematical Sunyaev-Zeldovich effects also by taking 
the motion of the Solar System into account. Numerical re- 
sults are presented in § 3, and finally, concluding remarks 
are given in § 4. 

2. Lorentz- boosted Kompaneets equation 

In the present section we extend the work of Nozawa, Itoh, 
& Kohyama (1998b) by first extending the Kompaneets 
equation to a system (the cluster of galaxies) moving 
with a peculiar velocity with respect to the CMBR. To 
formulate the kinetic equation for the photon distribu- 
tion function we use a relativistically covariant formalism 
(Berestetskii, Lifshitz, & Pitaevskii 1982; Buchler & Yueh 
1976). As a reference system, we have chosen the system 
which is fixed to the cosmic microwave background radia- 
tion (CMBR). The z-axis is fixed to a line connecting the 
observer and the center of mass of the cluster of galax- 
ies (CG), first assuming that the observer is fixed to the 
CMBR frame. We fix the positive direction of the z-axis 
as the direction of the propagation of a photon from the 
cluster to the observer. In this reference system, the cen- 
ter of mass of the CG is moving with a peculiar velocity 
/3(= v/c) with respect to the CMBR. For simplicity, we 
choose the direction of the velocity in the x-z plane, i.e. 

/3=(/3 x ,0,&). 

In the CMBR frame the four-momenta of the initial 
electron and photon arc p = (E,p) and k = (w,0,0, k), 
respectively, while the four-momenta of the final electron 
and photon arc p' = (E',p') and k' = (uj',k'), respec- 
tively. The electron distribution functions in the initial 
and final states are Fermi-like in the CG frame. They are 
related to the electron distribution functions in the CMBR 
frame as follows (Landau & Lifshitz 1975): 

f(E) = fc(E c ), (1) 

f(E') = f c (E' c ), (2) 

E c = 7 (E-p-p) , (3) 

E'c = 7 (E'-(3-p') , (4) 
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where the suffix C denotes the CG frame. 

In addition to 8 e , there is another parameter (3. For 
most of the CG, (3 <C 1 is realized. For example, (3 w 1/300 
for a typical value of the peculiar velocity u=l,000km/s. 
Therefore it should be sufficient to expand in powers of (3 
and to retain up to O{0 2 ) contributions. We assume the 
initial photon distribution of the CMBR to be Planckian 
with a temperature T : 

n (X) = ^—, (6) 



1 



where 
X 



k B T 



(7) 



Nozawa, Itoh & Kohyama (1998b) have obtained the 
following expression for the fractional distortion of the 
photon spectrum: 
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where # 7 is the angle between the directions of the pecu- 
liar velocity of the cluster (J3) and the initial photon mo- 
mentum (fc), which is chosen as the positive ^-direction. 
The reader should notice that this sign convention for the 
positive ^-direction is the opposite of the ordinary one. 
Thus a cluster moving away from the observer has f3 z < 0. 
However, because of the positive sign in front of the Pi((3 z ) 
term, one obtains An(X) < in this case, as one should. 
The coefficients are defined as follows: 
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X 

s 



Xcoth 

X 



sinh 



X 



(22) 
(23) 



where N e is the electron number density in the CG frame, 
the integral in Eq. (21) is over the photon path length in 
the cluster, and <tt is the Thomson cross section. 

In the next step of the calculation, we suppose that 
the observer's system (the Solar System) is moving with 
a velocity 0s(= vs/c) with respect to the CMBR. The 
photon distribution function in the CMBR system n(w, k) 
and the photon distribution function in the Solar System 
ns(tos,ks) are related by (Landau & Lifshitz 1975) 



n(to,k) = n s {us,k s ) , 

w = isi^s + 0s • k s ) ■ 
1 

7s = 



(24) 
(25) 

(26) 



The photon unit wave vector in the Solar System ks and 
the photon unit wave vector in the CMBR system k are 
related by (M0ller 1962) 
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where 0s is a un it vector in the direction of /3s- 

Therefore, our task is to rewrite Eq. (8) for the ex- 
pression corresponding to the Solar System. The isotropic 
Planck distribution in the CMBR system, Eq. (6), corre- 
sponds to the distribution including the dipolar distortion 
in the Solar System 
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9s being the angle between /3s and ks- Here we have used 
the relationship 



X = j s X s (l+0s^s)- 



(31) 



The velocity vecor j3 of the galaxy cluster in the CMBR 
system is transformed to the velocity vector f3' of the 
galaxy cluster in the Solar System by the relationship 
(M0ller 1962) 
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Therefore, Z in the CMBR system in Eq. (9) is trans- 
formed to 
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We write the photon distribution function in the Solar 
System frame that has a distortion caused by the thermal 
and kinematical Sunyaev-Zeldovich effects as 



n s (X s ,k s ) = n s {Xs,k s ) + An s (X Sl k s ). 



(34) 



Here Ans(Xs, ks) is given by 
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where one should rewrite X by using Eq. (31) and also 
0Pi0 z ) and P 2 P 2 {Pz) by using Eqs. (32), (33). The 
present result can be used in order to analyze the observa- 
tional data of the Sunyaev-Zeldovich effect. We note that 
0S = 1.241 x 10~ 3 is a known quantity and ^s~ cos9s is 
also a known quantity for the galaxy cluster under con- 
sideration. For present and near-future observations of the 
Sunyaev-Zeldovich effect, 2 terms in Eq. (35) are beyond 
observational accuracy, such that we only consider the first 
and third terms in the braces of Eq. (35). When the opti- 
cal depth of the galaxy cluster r is known, we obtain the 
best-fit values of 9 e and Z by fitting the observational 
data with the first and third terms of Eq. (35) written for 
the Solar System frame. In this way we can obtain the 
cluster temperature, as well as the cluster velocity in the 
direction of the photon propagation vector observed by 
the observer fixed to the CMBR system. 

In the somewhat more distant future accuracy in the 
observation of the Sunyaev-Zeldovich effect will be hope- 
fully improved dramatically. In that case one should also 
include the 2 terms in Eq. (35) and use the full expres- 
sion for analysis of the observational data. In doing so it 
should be possible to also determine the cluster velocity 
perpendicular to the photon propagation vector. We de- 
fine the distortion of the spectral intensity in the Solar 
System frame as follows: 



AI S = An s (Xs,ks)X 3 s . 



3. Numerical results 



(36) 



Equations (31) and (35) present our main results for the 
effect of the observer's the motion (the Solar System) . It 
should be emphasized here that because no expansions 
have been made so far in deriving Eqs. (31) and (35) for 
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Fig. 1. Spectral intensity of the thermal Sunyavev- 
Zeldovich effect for fcsT e = lOkeV and [is = —1 as a 
function of Xs- The solid curve is the full contribution. 
The dash-dotted curve is the lowest-order (Y ) contribu- 
tion. The dashed curve is the extraction of the full (3s 
corrections (multiplied by 100). The dotted curve shows 
results of Chluba ct al. (2005) (multiplied by 100). Refer 
to the main text for detailed explanations. 

Psi they are general expressions for arbitrary values of 
(3 s- In practical cases, however, (3 s -C 1 is achieved. By 
expanding X around Xs and keeping only lowest order 
terms for both (3s and 9 e , we have reproduced the Eqs. 
(12a) and (12b) of Chluba, Huetsi, and Sunyacv (2005). 

In Fig. 1, we plot the spectral intensity of the thermal 
Sunyavev-Zeldovich effect for fcgT e = lOkeV and us = — 1 
as a function of Xs- Here the solid curve is the full con- 
tribution (Y"o to Y4 terms) of the first line of Eq. (35), and 
the dash-dotted curve is the lowest-order (Yo term) con- 
tribution of the first line of Eq. (35). The dashed curve 
is the extraction of the full (3s corrections of the first line 
of Eq. (35), which is multiplied by 100 in order to be 
visible in the same figure. This curve contains full (3s cor- 
rections. The dotted curve shows the results of Chluba, 
Huetsi, and Sunyaev (2005) (Eq. (12a) in their paper), 
where only lowest-order corrections are included for both 
(3s and 9 e . This correction is also multiplied by 100 in 
order to be visible in the same figure. It should be noted 
that the (3 S corrections are totally negligible. The essential 
difference in two curves are due to higher-order relativistic 
corrections of 9 ei which is present in Eq. (35) but is ab- 
sent in Eq. (12a) of Chluba, Huetsi, and Sunyaev (2005). 
Therefore, higher order relativistic corrections of 9 e are 
important in discussing the effect of the observer's motion 
in the thermal Sunyaev-Zeldovich effect with high preci- 
sion. 

In Fig. 2, we have also plotted the spectral intensity 
of the kincmatical Sunyavev-Zeldovich effect for ksT e = 
lOkeV, (3 = (3 Z — 1/300, and /is = -1 as a function of Xs- 
Note that the velocity of the cluster relative to the CMBR 
and the velocity of the observer relative to the CMBR are 
parallel in this case. The solid curve is the full contribu- 



kinematical SZ 10keV 



/3=/? z = 1/300 




/? s corrections (x100) 



5 10 15 

Fig. 2. Spectral intensity of the kincmatical Sunyavev- 
Zeldovich effect for k B T e = lOkeV, (3 = (3 Z = 1/300, and 
/is = —1 as a function of Xs- The solid curve is the full 
contribution. The dash-dotted curve is the lowest-order 
contribution (3P\{(3 Z ). The dashed curve is the extraction 
of the full (3s corrections (multiplied by 100). The dotted 
curve shows the results of Chluba et al. (2005) (multiplied 
by 100). Refer to the main text for a detailed explanations. 

tion (1 to C2 terms) of the third line of Eq. (35), and the 
dash-dotted curve is the lowest-order contribution of the 
third line of Eq. (35). The dashed curve is the extraction of 
the full (3s corrections of the third line of Eq. (35), which 
is multiplied by 100 in order to be visible in the same fig- 
ure. This curve contains full (3s corrections. The dotted 
curve shows the results of Chluba, Huetsi, and Sunyaev 
(2005) (Eq. (12b) in their paper), where only lowest-order 
corrections are included for both (3s and 9 e . This correc- 
tion is also multiplied by 100 in order to be visible in the 
same figure. It should also be noted that the (3 S correc- 
tions are totally negligible. The essential difference in two 
curves are due to higher-order relativistic corrections of 
9 e , which is present in Eq. (35) but is absent in Eq. (12b) 
of Chluba, Huetsi, and Sunyaev (2005). Therefore, higher 
order relativistic corrections of 9 e are again important in 
discussing the effect of the motion of the observer in the 
kinematical Sunaev-Zeldovich effect with high precision. 

4. Concluding remarks 

We have calculated the effect of the motion of the ob- 
server (the Solar System) on the relativistic thermal and 
kinematical Sunyaev-Zeldovich effects with the power se- 
ries approximation in terms of 9 e = kBT e /mc 2 . We made 
full use of the Lorentz covariance properties of the prob- 
lem to obtain solutions with high precision, and con- 
firmed the correctness of the results recently obtained 
by Chluba, Huetsi, and Sunyaev in the lowest-order of 
the observer's velocity (3s- We have given a more gen- 
eral analytic expression for the thermal and kinematical 
Sunyaev-Zeldovich effects corresponding to the observer's 
system (the Solar System). We also found that the ef- 
fect of the motion of the observer (the Solar System) on 
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the thermal and kincmatical Sunyaev-Zeldovich effects is 
marginally important for present and near-future observa- 
tions. However, for the next-generation high-precision ob- 
servations, we find the effect of the motion of the observer 
(the Solar System) could be important and should be 
taken into account for analysis of the observational data. 
We also note that the next-generation high-precision kine- 
matic Sunyaev-Zeldovich effect observations will hopefully 
allow us to determine the cluster velocity perpendicular to 
the photon propagation vector, as well as the cluster ve- 
locity in the direction of the photon propagation vector. 
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